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Abstract: 


This  research  report  includes  three  parts.  First,  preparation  of  a  high-concentration  Au  nanoring  (NR) 
water  solution  and  its  applications  to  the  enhancement  of  image  contrast  in  optical  coherence 
tomography  (OCT)  and  the  generation  of  photothermal  effect  in  a  sample  through  localized  surface 
plasmon  (LSP)  resonance  are  demonstrated.  Au  NRs  are  first  fabricated  on  a  sapphire  substrate  with 
colloidal  lithography  and  secondary  sputtering  of  Au,  and  then  transferred  into  water  solution  through 
a  liftoff  process.  By  controlling  the  NR  geometry,  the  LSP  dipole  resonance  wavelength  in  tissue  can 
cover  the  spectral  range  of  1300  nm  for  OCT  scanning  of  deep  tissue  penetration.  The  extinction  cross 
sections  of  the  fabricated  Au  NRs  in  water  are  estimated  to  give  the  levels  of  10'10-10'9  cm2  near  their 
LSP  resonance  wavelengths.  The  fabricated  Au  NRs  are  then  delivered  into  pig  adipose  samples  for 
OCT  scanning.  It  is  observed  that  when  resonant  Au  NRs  are  delivered  into  such  a  sample,  LSP 
resonance-induced  Au  NR  absorption  results  in  a  photothermal  effect,  making  the  opaque  pig  adipose 
cells  transparent.  Also,  the  delivered  Au  NRs  in  the  intercellular  substance  enhance  the  image  contrast 
of  OCT  scanning  through  LSP  resonance-enhanced  scattering.  By  continuously  OCT  scanning  a 
sample,  both  photothermal  and  image  contrast  enhancement  effects  are  observed.  However,  by 
continually  scanning  a  sample  with  a  low  scan  frequency,  only  the  image  contrast  enhancement  effect 
is  observed. 

Second,  the  simulation  results  of  absorption  enhancement  in  an  amorphous-Si  (a-Si)  solar  cell  by 
depositing  metal  nanoparticles  (NPs)  on  the  device  top  and  embedding  metal  NPs  in  a  layer  above  the 
A1  back-reflector  are  demonstrated.  The  absorption  increase  results  from  the  near-field  constructive 
interference  of  electromagnetic  waves  in  the  forward  direction  such  that  an  increased  amount  of 
sunlight  energy  is  distributed  in  the  a-Si  absorption  layer.  Among  the  three  used  metals  of  Al,  Ag,  and 
Au,  Al  NPs  show  the  most  efficient  absorption  enhancement.  Between  the  two  used  NP  geometries, 
Al  nanocylinder  (NC)  are  more  effective  in  absorption  enhancement  than  Al  nanosphere  (NS).  Also,  a 
random  distribution  of  isolated  metal  NCs  can  lead  to  higher  absorption  enhancement,  when  compared 
with  the  cases  of  periodical  metal  NC  distributions.  Meanwhile,  the  fabrication  of  both  top  and 
bottom  Al  NCs  in  a  solar  cell  results  in  further  absorption  enhancement.  Misalignments  between  the 
top  and  bottom  Al  NCs  do  not  significantly  reduce  the  enhancement  percentage.  With  a  structure  of 
vertically  aligned  top  and  bottom  Al  NCs,  solar  cell  absorption  can  be  increased  by  52  %. 

Third,  the  simulation  results  of  the  enhanced  radiated  power  of  a  radiating  dipole  when  it  is 
coupled  with  the  localized  surface  plasmon  resonance  (LSPR)  of  a  nearby  metal  nanosphere  (NS)  by 
including  the  influence  of  the  LSPR  on  the  strength  of  the  dipole,  which  is  represented  by  a  two-level 
system  in  our  semi-classical  model,  are  demonstrated.  It  is  shown  that  with  this  feedback  influence, 
the  total  radiated  power  can  be  further  enhanced  by  more  than  several  tens  percents,  when  compared 
with  the  condition  of  ignoring  the  feedback  effect.  In  comparing  the  results  of  two  different  dipole 
orientations,  it  is  discovered  that  at  a  fixed  distance  from  the  metal  NS,  the  radial  dipole  orientation 
leads  to  more  significant  enhancement  of  total  radiated  power,  when  compared  to  the  orbital  dipole 
orientation. 
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Part  I:  Au  Nanorings  for  Enhancing  Absorption  and  Backscattering  at  Localized  Surface 
Plasmon  Resonance 

1.  Introduction 

With  its  unique  scattering  and  absorption  features  near  resonance,  localized  surface  plasmon  (LSP) 
induced  on  a  metal  nanoparticle  (NP)  has  found  broad  applications  in  emission  and  absorption 
enhancements  [1-5],  biomedical  imaging  contrast  improvement  [6-10],  and  photothermal  therapy 
[11-14].  Among  various  metals  for  fabricating  NPs  for  biomedical  applications,  Au  is  preferred  due  to 
its  biocompatibility.  By  connecting  the  anti-body  of  targeted  cells  to  Au  NPs  through  a 
bio-conjugation  process,  the  NPs  can  attach  the  targeted  cells  and  possibly  be  internalized  into  cells 
when  they  are  delivered  into  bio-tissue.  In  optical  imaging  of  bio-tissue,  the  strong  scattering  of  the 
Au  NPs  around  the  LSP  resonance  wavelength  can  lead  to  a  higher  image  contrast  of  the  targeted  cells. 
Also,  the  strong  LSP  absorption  results  in  local  heating  in  the  tissue  and  can  cause  the  death  of 
targeted  cells,  leading  to  the  effect  of  photothermal  therapy.  For  such  biomedical  applications,  the 
widely  used  Au  NPs  include  nanorods  [15-17],  nanoshells  [18-20],  and  nanocages  [21-23].  However, 
most  of  the  LSP  resonance  wavelengths  of  those  reported  Au  NPs  are  shorter  than  1000  nm.  From  the 
viewpoint  of  optical  imaging  application,  particularly  of  optical  coherence  tomography  (OCT), 
although  a  shorter  wavelength  around  800  nm  can  lead  to  higher  imaging  resolution,  the  use  of  a  light 
source  of  a  longer  wavelength  around  1300  nm  can  result  in  deeper  tissue  penetration  due  to  weaker 
tissue  scattering.  The  issue  of  deeper  tissue  penetration  of  light  is  important  not  only  for  imaging  but 
also  for  photothermal  therapy.  Although  an  Au  nanorod  of  a  larger  aspect  ratio  can  push  the  LSP 
resonance  wavelength  into  the  spectral  range  around  1300  nm,  it  is  usually  difficult  to  control  the 
uniformity  of  nanorods  and  hence  the  application  becomes  unreliable.  Also,  although  theoretically 
nanoshells  and  nanocages  of  thin  Au  layers  can  lead  to  LSP  resonances  at  longer  wavelengths,  the 
non-uniform  shell/wall  thickness  results  in  significantly  broadened  and  hence  weaker  resonance  peaks. 
Reliable  Au  NPs  of  other  geometries  for  longer-wavelength  and  stronger  LSP  resonance  are  needed 
for  further  developments  of  targeted  optical  imaging  and  photothermal  therapy.  For  this  purpose,  Au 
nanoring  (NR)  shows  the  advantage  of  long-wavelength  LSP  resonance.  The  fabrication  of  Au  NRs  on 
a  substrate  has  been  reported  based  on  the  methods  of  colloidal  lithography  and  secondary  sputtering 
of  metal  [24-26].  However,  the  delivery  of  such  Au  NRs  into  tissues  or  cells  for  optical  imaging, 
therapy,  and  other  biomedical  applications  has  not  yet  been  reported. 

In  this  part  of  report,  we  report  the  fabrication  of  Au  NRs  on  sapphire  substrate  using  the 
methods  of  colloidal  lithography  and  secondary  sputtering  of  Au  and  the  development  of  a  process  for 
NR  liftoff  from  the  substrate.  The  aqueous  solution  of  Au  NRs  is  prepared  for  measuring  the  LSP 
resonance  properties.  The  results  agree  reasonably  well  with  those  of  numerical  simulations.  The  Au 
NRs  are  then  delivered  into  pig  adipose  samples  for  OCT  scanning  to  demonstrate  the  enhanced 
absorption  and  scattering  behaviors  of  LSP  resonance.  It  is  observed  that  when  the  LSP  resonance 
wavelength  coincides  with  the  OCT  light  source  spectrum,  the  resonant  Au  NRs  lead  to  the  heating  of 
adipose  cells  and  higher  image  contrast.  Such  features  are  not  observed  when  another  Au  NR  sample, 
in  which  the  LSP  resonance  wavelength  significantly  deviates  from  the  OCT  light  source  spectrum,  is 
delivered  into  a  pig  adipose  sample.  OCT  is  a  widely  used  optical  imaging  technique  for  the  diagnoses 
of  many  diseases  [27-29],  including  early-diagnosis  of  oral  cancer  and  precancer  [30-32].  It  is  based 
on  the  principle  of  low-coherence  interferometry.  In  most  human  tissue,  OCT  scanning  can  lead  to 
high-resolution  imaging  with  a  depth  of  2-3  mm,  showing  its  most  advantageous  feature  when 
compared  with  optical  microscopy  techniques.  With  a  light  source  of  >100  nm  in  spectral  width 
within  the  spectral  window  of  700-1350  nm,  the  tissue  resolution  can  be  higher  than  7  microns.  With 
the  interferometry  nature  of  OCT  and  the  coherent  property  of  LSP-induced  scattering  of  Au  NPs, 
OCT  represents  an  effective  approach  for  monitoring  the  LSP-induced  functions  in  bio-tissue. 

This  part  of  report  is  organized  as  follows:  In  Section  2,  the  procedures  of  fabricating  Au  NRs  on 
sapphire  substrate  and  then  in  aqueous  solution  are  described.  The  measured  and  simulated  LSP 
resonance  properties  of  the  fabricated  Au  NRs  are  demonstrated  in  Section  3.  Next,  OCT  scanning 
images  created  by  delivering  Au  NRs  into  pig  adipose  samples  are  shown  and  discussed  in  Section  4. 
Finally,  conclusions  are  drawn  in  Section  5. 

2.  Fabrication  of  Au  Nanorings 

The  fabrication  of  Au  NRs  starts  with  a  colloidal  lithography  process  through  electrostatic  adsorption 
of  polystyrene  colloidal  particles  (Duke  Scientific  Corporation)  on  a  Si02  layer  of  20  nm  in  thickness 
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[33],  which  is  coated  on  a  sapphire  substrate  with  plasma-enhanced  chemical  vapor  deposition.  Two 
samples  are  prepared  with  colloidal  particle  diameters  at  150  (sample  A)  and  100  (sample  B)  nm. 
Then,  an  Au  film  -30  nm  in  thickness  is  deposited  on  top  of  the  substrate  to  serve  as  the  source  of 
secondary  sputtering.  Next,  reactive  ion  etching  (RIE)  is  applied  by  using  CHF3  as  etchant  for  forming 
Au  ring  structures  around  the  colloidal  particles  through  the  processes  of  secondary  sputtering  and 
removing  the  rest  Au  atoms.  During  the  RIE  process,  the  colloidal  particles  are  partially  etched.  An  02 
plasma  etching  procedure  completely  cleans  the  colloidal  particle  material,  leaving  Au  NRs  on 
Si02-deposited  sapphire  substrate.  A  typical  scanning  electron  microscopy  (SEM)  image  of  the  Au 
NRs  of  sample  A  on  the  Si02-deposited  sapphire  substrate  is  shown  in  Fig.  1.  Here,  the  average  outer 
diameter  of  those  Au  NRs  is  about  174  nm  with  the  inner  diameter  at  —150  nm  since  Au  atoms  are 
sputtered  onto  the  sidewall  of  a  colloidal  particle.  Therefore,  the  average  thickness  of  the  NRs  is  about 
12  nm.  The  average  outer  diameter  of  the  Au  NRs  of  sample  B  is  about  118  nm  with  the  average  Au 
thickness  at  -9  nm. 

The  Au  NRs  are  then  surface  modified  by  APTS  ((3-aminopropyl)  trimethoxysilane,  97%,  Alfa 
Aesar)  to  prevent  them  from  aggregation  in  water.  In  this  process,  a  1:1  water-ethanol  solution  (100 
mL)  is  combined  with  50  pL  of  APTS,  and  is  maintained  at  pH  10  by  adding  NH3  solution.  The 
substrate  samples  are  immersed  in  this  mixture  for  more  than  24  hours  and,  subsequently,  are  rinsed 
with  ethanol  and  dried  with  N2  flow.  To  transfer  the  substrate-attached  Au  NRs  into  water  solution, 
the  substrate  samples  are  immersed  in  a  KOH  solution  (85%,  13.6%  w.t.  in  water)  for  15  min  to  etch 
the  Si02  layer  between  Au  NRs  and  sapphire  substrate,  followed  by  a  step  of  gentle  rinse  with 
de-ionized  water  to  remove  residual  KOH.  Then,  the  substrate  sample  is  placed  in  a  glass  bottle  filled 
with  the  proper  amount  of  de-ionized  water.  After  sonication,  the  Au  NRs  are  transferred  into  water 
solution.  To  increase  the  stability  of  NRs,  PEG-Thiol  (mPEG-Thiol-5000,  Laysan  Bio  Inc.,  Arab,  AL) 
is  added  to  the  NR  solution  to  reach  a  concentration  of  2.5 xlCT5  M.  The  mixture  is  then  incubated  at 
room  temperature  for  20  hours.  To  deliver  the  NRs  into  tissue  for  OCT  scanning,  the  NR  solution  is 
centrifuged  to  give  a  higher  particle  concentration.  Figure  2  shows  a  plan- view  SEM  image  of  Au  NRs 
of  sample  A  after  liftoff.  The  SEM  image  was  taken  when  the  dried  Au  NRs  are  placed  on  a  Si 
substrate.  Here,  one  can  see  that  the  inner  and  outer  diameters  of  the  NRs  change  slightly  along  their 
ring  axes. 

3.  Surface  Plasmon  Characteristics  of  Au  Nanorings 

Optical  transmissions  of  the  water  solutions  of  Au  NRs  are  measured  with  a  UV-vis 
spectrophotometer  (Jasco  V-570)  to  show  the  normalized  extinction  spectra  in  Fig.  3  (the  left  ordinate). 
Here,  one  can  see  the  major  peaks  of  extinction,  corresponding  to  the  LSP  dipole  resonances  in  water, 
at  1240  and  1030  nm  for  samples  A  and  B,  respectively.  In  this  figure,  the  spectrum  of  the 
sweeping-frequency  laser  source  of  the  used  OCT  system  is  also  shown  for  comparison  (the  right 
ordinate).  The  light  source  spectrum  with  its  peak  at  1310  nm  and  spectral  full-width  at  half-maximum 
at  -100  nm  is  covered  by  the  long-wavelength  side  of  the  LSP  dipole  resonance  range  of  sample  A. 
However,  the  LSP  dipole  resonance  range  of  sample  B  is  quite  far  away  from  the  OCT  light  source 
spectrum.  The  secondary  peaks  of  the  extinction  spectra  in  Fig.  3  correspond  to  the  mixture  of  the 
higher-order  resonance  across  the  ring  and  the  resonance  along  the  ring  axis  (see  the  discussions  for 
Figs.  4  and  5). 

To  further  understand  the  resonance  behaviors  of  the  fabricated  Au  NRs,  numerical  simulations 
based  on  the  finite- element  method  are  performed  by  assuming  an  ideal  ring  shape  shown  in  the  insert 
of  Fig.  4.  In  Fig.  4,  the  extinction  spectra  of  the  Au  NRs  are  obtained  with  the  excitation  light  incident 
along  the  z  axis  and  incident  linear  polarization  in  the  y  direction.  Here,  the  normalized  extinction 
spectra  of  the  two  samples  with  ambient  refractive  indices  at  1.33  (corresponding  to  NRs  in  water  - 
curves  A(1.33)  and  B(1.33))  and  1.4  (corresponding  to  NRs  in  tissue  -  curves  A(1.4)  and  B(1.4))  are 
shown.  The  used  NR  parameters  include  the  outer  diameter  at  179  (120)  nm,  inner  diameter  at  149 
(97)  nm,  and  height  at  95  (60)  nm  for  sample  A  (B).  For  simulations,  the  dispersive  dielectric  constant 
of  Au  follows  the  experiment  data  in  ref.  [34].  Although  the  SEM  image  in  Fig.  2  shows  that  the 
diameter  of  an  NR  changes  along  its  axis,  the  numerical  results  in  Fig.  4  show  the  reasonable 
agreements  of  LSP  resonance  spectra  with  the  measurements.  In  Fig.  4,  one  can  also  see  that  by 
increasing  the  ambient  refractive  index,  the  LSP  resonance  wavelengths  are  red  shifted.  In  this 
situation,  the  spectral  range  of  the  LSP  dipole  resonance  of  sample  A  completely  covers  the  OCT  light 
source  spectrum.  Meanwhile,  that  of  sample  B  only  slightly  overlaps  with  the  OCT  light  source 
spectrum.  Therefore,  it  is  expected  that  only  sample  A  can  significantly  induce  LSP 
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resonance-enhanced  backscattering  and  absorption  in  bio-tissue.  It  is  noted  that  the  extinction  spectra 
obtained  from  water  solution  measurements  shown  in  Fig.  3  are  significantly  broader  than  their 
counterparts  obtained  from  single-NR  simulations  shown  in  Fig.  4.  This  difference  can  be  attributed  to 
(1)  the  non-uniform  shape  and  size  distributions  among  different  NRs,  (2)  the  random  orientation 
distribution  of  Au  NRs  in  water  solution,  and  (3)  the  possible  LSP  coupling  between  neighboring  NRs 
in  water  solution. 

For  further  understanding  the  LSP  resonance  properties  of  the  fabricated  Au  NRs  in  water 
solution,  we  also  evaluate  the  extinction  spectra  of  LSP  resonances  across  the  ring  and  along  the  ring 
axis  by  assuming  that  the  excitation  light  is  incident  along  the  y-axis  and  the  incident  linear 
polarizations  are  in  the  x  and  z  directions,  respectively.  Figure  5  shows  the  simulation  results  of  the 
two  samples,  including  the  resonances  across  the  ring  (curves  A/ring  and  B/ring)  and  those  along  the 
ring  axis  (curves  A/axis  and  B/axis).  In  the  curves  for  the  resonances  across  the  ring,  one  can  see  the 
major  peaks  of  dipole  resonance  the  same  as  those  shown  in  Fig.  4.  The  minor  peaks  at  shorter 
wavelengths  correspond  to  higher-order  resonance  modes.  The  combination  of  the  higher-order 
cross-ring  resonance  mode  with  LSP  modes  along  the  ring  axis  corresponds  to  the  secondary  peaks 
between  500  and  900  nm  in  Fig.  3.  The  data  shown  in  Fig.  3  represent  the  total  resonance  effect  of 
numerous  Au  NRs  of  random  orientation  distributions  in  water. 

The  accurate  evaluation  of  NP  concentration  in  solution  still  represents  a  major  problem  in  the 
development  of  Au  NP  application.  The  NR  concentration  in  water  solution  can  be  estimated  by  first 
counting  the  Au  NR  surface  density  on  sapphire  substrate  based  on  the  statistics  of  SEM  observations 
to  give  1.44  x  109  cm'2  for  sample  A  and  2.64  x  109  cm'2  for  sample  B.  This  surface  density  is 
multiplied  by  the  substrate  area  (a  quarter  of  a  2-inch  sapphire  substrate)  to  give  the  total  NR  numbers 
on  a  substrate  at  7.29  x  109  NRs  for  sample  A  and  1.34  x  1010  NRs  for  sample  B.  These  Au  NRs  are 
transferred  into  a  water  solution  of  4  mL  in  volume  to  give  particle  concentrations  of  1.82  x  109  and 
3.35  x  109  cm'3  for  samples  A  and  B,  respectively.  It  is  noted  that  in  this  concentration  estimation 
procedure,  errors  can  be  made  through  inaccuracy  in  evaluating  the  NR  surface  density  on  substrate 
and  NR  loss  in  transferring  NRs  into  water  solution.  Therefore,  the  NR  concentrations  of  samples  A 
and  B  are  estimated  to  be  -1.09  x  109  and  -2.01  x  109  cm'3,  respectively,  by  assuming  that  40  %  NRs 
are  lost  in  the  process  of  transferring  NRs  into  water  solution.  The  water  solutions  of  the  two  samples 
are  then  condensed  by  a  factor  of  about  12  to  give  the  final  concentrations  of -1.31  x  1010  and  -2.41  x 
1010  cm'3  in  samples  A  and  B,  respectively,  for  OCT  scanning  application.  We  also  use  the  method  of 
inductively  coupled  plasma  mass  spectrometry  for  estimating  the  NR  concentrations  to  give  1.07  x 
1010  and  3.58  x  1010  cm'3  for  samples  A  and  B,  respectively,  after  condensation.  Considering  all  the 
inaccuracies  in  measurements  and  estimations,  we  suppose  that  these  results  agree  reasonably  well 
with  the  aforementioned  values  obtained  from  surface  density  estimation. 

The  extinction  cross  sections  of  the  NRs  can  be  evaluated  based  on  transmission  measurement. 
By  using  the  NR  solutions  before  condensation,  we  obtained  the  transmittances  of  54.58  and  56.47  % 
through  the  solutions  of  1  cm  in  length  for  samples  A  and  B,  respectively,  at  their  individual  resonance 
wavelengths,  i.e.,  at  1240  and  1030  nm  for  samples  A  and  B,  respectively.  From  the  measurement 
results,  we  can  obtain  the  NR  extinction  cross  sections  of  5.56  x  10'10  and  2.84  x  10'10  cm2  for  samples 
A  and  B,  respectively,  at  their  individual  resonances.  Figure  6  shows  the  extinction  (Ext),  scattering 
(Sea),  and  absorption  (Abs)  cross  sections  of  single  Au  NRs  in  water  as  functions  of  wavelength  for 
samples  A  and  B  obtained  from  numerical  simulations.  Here,  one  can  read  the  extinction  cross  sections 
of  3.13  x  10'9  and  1.70  x  10'9  cm2  at  individual  resonance  wavelengths  for  samples  A  and  B, 
respectively.  The  simulation  results  of  single  NRs  are  about  5.8  times  those  from  transmission 
measurement.  Several  reasons  can  be  used  for  explaining  this  difference,  including  (1)  the 
non-uniform  size  and  geometry  distributions  of  NRs  (leading  to  the  reduction  of  overall  resonance 
strength  and  broadening  of  resonance  spectrum),  (2)  the  random  orientation  distribution  of  NRs 
(resulting  in  weaker  overall  resonance  strength),  (3)  the  difficulty  of  accurately  evaluating  NR 
concentration  (underestimating  the  extinction  cross  section  based  on  transmission  measurement),  and 
(4)  the  possibility  of  LSP  coupling  between  neighboring  NRs  (also  leading  to  the  reduction  of  overall 
resonance  strength  and  broadening  of  resonance  spectrum).  In  Fig.  6,  it  is  noted  that  the  spectral  peaks 
of  scattering  and  absorption  cross  sections  are  slightly  different  in  both  samples.  Also,  scattering  is 
significantly  stronger  than  absorption  in  sample  A.  However,  scattering  and  absorption  strengths 
become  comparable  in  sample  B. 

4.  Optical  Coherence  Tomography  Images 
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In  the  used  OCT  system,  a  sweeping-frequency  laser  source  with  the  output  spectrum  shown  in  Fig.  3 
is  used  as  the  light  source.  The  light  source  is  connected  to  a  Mech-Zehnder  interferometer,  which 
consists  of  two  optical  couplers  and  two  circulators.  The  interference  fringe  signal  is  detected  by  a 
balanced  photo-detector  and  sampled  by  a  high-speed  digitizer.  The  achieved  system  sensitivity  and 
axial  resolution  in  tissue  are  103  dB  and  6  m  at  the  depth  of  1  mm,  respectively.  The  detailed 
descriptions  about  the  OCT  system  and  its  operation  condition  can  be  found  in  ref.  [32]. 

Although  the  photothermal  effect  generated  by  LSP  resonance  of  Au  NP  has  been  widely 
observed,  the  measurement  of  local  temperature  variation  in  a  bio-sample  is  still  difficult.  To 
demonstrate  the  effects  of  LSP  resonance-enhanced  scattering  and  absorption  in  bio-tissue,  we  chose 
pig  adipose  as  our  sample  for  OCT  scanning.  Pig  adipose  cells  of  100-200  m  in  size  are  opaque  in 
the  visible-near-infrared  range  at  a  temperature  lower  than  25  °C.  When  such  a  cell  is  heated  to  reach  a 
higher  temperature,  it  becomes  transparent  [35].  Therefore,  pig  adipose  can  be  used  for  demonstrating 
the  local  heating  effect  through  LSP  resonance.  In  our  experiment,  fresh  pig  adipose  samples  are 
refrigerated  before  OCT  scanning.  After  removing  a  sample  from  the  refrigerator,  it  takes  about  four 
hours  to  observe  transparent  cells  at  the  ambient  temperature  of  25  °C.  Figures  7(a)-7(d)  show  the 
OCT  scanning  images  at  the  same  location  of  a  pig  adipose  sample  before,  0  min,  30  min,  and  60  min 
after  the  application  of  a  water  solution  droplet  of  non-resonant  Au  NRs  (sample  B)  onto  the  sample 
surface,  respectively.  The  OCT  scanning  experiment  was  started  right  after  the  sample  was  moved  out 
from  the  refrigerator.  The  sample  was  continuously  scanned  with  the  OCT  system  (about  15  scans  per 
sec).  In  Fig.  7(b),  one  can  clearly  see  the  backscattering  of  Au  NRs  in  the  water  solution  droplet.  The 
broad  resonance  spectrum  of  sample  B  shown  in  Fig.  3  partially  overlaps  with  the  OCT  light  source 
spectrum,  leading  to  observable  LSP  resonance-induced  backscattering  in  OCT  scanning.  After  the 
non-resonant  Au  NRs  are  delivered  into  the  adipose  sample,  as  shown  in  Figs.  7(c)  and  7(d),  one  can 
hardly  observe  any  change  in  the  OCT  image  from  that  before  the  application  of  non-resonant  Au 
NRs. 

The  OCT  scanning  results  with  resonant  Au  NRs  are  quite  different.  Figures  8(a)-8(d)  show  the 
OCT  scanning  images  at  the  same  location  of  a  pig  adipose  sample  before,  0  min,  30  min,  and  60  min 
after  the  application  of  a  water  solution  droplet  of  resonant  Au  NRs  (sample  A)  onto  the  sample 
surface,  respectively.  The  sample  was  continuously  scanned  with  the  OCT  system.  In  Fig.  8(b),  one 
can  see  the  significantly  stronger  backscattering  signal  from  the  water  solution  droplet,  when 
compared  with  that  in  Fig.  7(b).  By  observing  Figs.  8(c)  and  8(d)  closely,  dark-hole  structures  of 
100-200  m  in  size  (transparent  adipose  cells)  can  be  seen.  The  adipose  images  in  Figs.  8(c)  and  8(d) 
are  unclear  because  the  Au  NRs  accumulated  on  the  sample  surface  block  the  incident  light  in  OCT 
scanning,  as  demonstrated  by  the  thick  bright  surface  boundaries  in  the  figures.  Figure  8(e)  shows  the 
OCT  image  around  60  min  after  the  application  of  Au  NR  solution  when  a  pure  water  droplet  is  used 
to  wash  out  the  accumulated  Au  NRs  on  the  surface.  Here,  clear  hole  structures  showing  transparent 
adipose  cells  can  be  seen,  indicating  the  local  heating  effect  due  to  LSP  resonance-induced  absorption 
of  Au  NRs.  Since  the  cells  look  so  transparent,  it  is  believed  that  most  of  the  Au  NRs  are  delivered 
into  the  intercellular  substance  for  enhancing  the  cell  imaging  contrast  in  OCT  scanning.  The  image 
contrast  enhancement  with  resonant  Au  NRs  can  be  observed  by  comparing  Figs.  9(a)  with  9(b). 
Figure  9(a)  is  a  duplicate  of  Fig.  7(e)  except  that  a  horizontal  dashed  line  is  drawn  for  OCT  signal 
intensity  demonstration  (see  Fig.  10).  The  OCT  scanning  image  in  Fig.  9(b)  is  obtained  by  scanning  a 
fresh  pig  adipose  sample  (without  being  refrigerated  or  using  any  Au  NR).  Here,  although  the  hole 
structures  can  be  clearly  seen  in  Fig.  9(b),  its  image  contrast  is  lower.  The  OCT  lateral  line-scan  signal 
intensity  distributions  along  the  horizontal  dashed  lines  in  Figs.  9(a)  and  9(b)  are  plotted  in  Fig.  10. 
Here,  one  can  clearly  see  the  enhanced  image  contrast  after  Au  NRs  are  delivered  into  the  adipose 
sample.  Therefore,  one  can  realize  the  effect  of  enhanced  scattering  due  to  LSP  resonance  besides  the 
photothermal  effect  for  making  the  adipose  cell  transparent. 

The  results  shown  in  Figs.  7-10  are  obtained  from  the  experiments  of  continuously  scanning  the 
samples  with  the  OCT  system.  In  each  OCT  scan,  the  incident  light  can  produce  heating  due  to 
LSP-resonance-induced  absorption  enhancement.  Therefore,  we  have  observed  a  combined  effect  of 
enhanced  absorption  and  scattering.  To  independently  observe  the  LSP-resonance-induced  scattering 
enhancement,  we  can  minimize  the  photothermal  effect  by  continually  scanning  the  sample  with  a  low 
scan  frequency.  Figures  11  (a)- 11(d)  show  the  OCT  images  of  a  pig  adipose  sample  taken  before,  0 
min,  30  min,  and  60  min  after  the  application  of  a  droplet  of  Au  NR  water  solution  of  sample  A, 
respectively,  when  the  sample  is  scanned  once  every  30  min  such  that  the  photothermal  effect  is 
minimized.  Figure  11(e)  shows  the  image  similar  to  that  in  Fig.  11(d)  after  the  Au  NRs  accumulated 
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on  the  surface  are  washed  out.  The  pig  adipose  sample  was  refrigerated  before  the  OCT  scanning 
experiment.  Here,  one  can  see  that  adipose  cells  remain  opaque  through  the  60-min  period  of 
experiment,  implying  that  the  photothermal  effect  is  essentially  deleted.  Although  the  adipose  cells  are 
opaque,  their  boundaries  can  still  be  seen  with  OCT  scanning.  Such  boundaries  correspond  to  the 
sharp  peaks  in  Fig.  12,  in  which  two  OCT  line-scan  intensity  distribution  curves  along  the  vertical 
lines  in  Figs.  11(a)  (Before)  and  11(b)  (Wash  Out)  are  shown.  Here,  one  can  clearly  see  that  after  the 
delivery  of  resonant  Au  NRs  into  the  intercellular  substance  of  the  adipose  sample,  the  image  contrast 
is  enhanced  under  the  condition  of  weak  photothermal  effect.  The  results  shown  in  Figs.  11  and  12 
clearly  demonstrate  the  effect  of  enhanced  image  contrast.  It  is  noted  that  the  intensity  levels  in  Figs. 
10  and  12  are  obtained  from  different  image  process  procedures  and  their  comparison  does  not  make 
any  sense. 

5.  Conclusions 

In  summary,  the  preparation  of  Au  NR  water  solution  of  high  concentration  and  their  applications  to 
the  enhancement  of  image  contrast  in  OCT  and  the  generation  of  photothermal  effect  in  bio-tissue 
have  been  demonstrated.  Au  NRs  were  first  fabricated  on  sapphire  substrate  with  colloidal  lithography 
and  secondary  sputtering  of  Au,  and  then  transferred  into  the  water  solution  through  a  liftoff  process. 
By  controlling  the  NR  geometry,  the  LSP  dipole  resonance  wavelength  in  tissue  could  cover  the 
spectral  range  around  1300  nm  for  OCT  scanning  of  deep  tissue  penetration.  The  extinction  cross 
sections  of  the  fabricated  Au  NRs  in  water  at  their  LSP  resonance  wavelengths  were  calibrated  to  give 
an  order  of  magnitude  of  10'10-10'9  cm2.  The  fabricated  Au  NRs  were  delivered  into  pig  adipose 
samples  for  OCT  scanning.  It  was  observed  that  when  resonant  Au  NRs  were  delivered  into  such  a 
sample,  LSP  resonance-induced  Au  NR  absorption  led  to  the  photothermal  effect  of  making  the 
opaque  pig  adipose  cells  transparent.  Also,  the  delivered  Au  NRs  in  the  intercellular  substance  could 
enhance  the  image  contrast  of  OCT  scanning  through  LSP  resonance-enhanced  scattering.  By 
continuously  OCT  scanning  an  adipose  sample,  both  photothermal  and  image  contrast  enhancement 
effects  could  be  observed.  However,  by  continually  OCT  scanning  a  sample  with  a  low  scan  frequency, 
only  the  image  contrast  enhancement  effect  could  be  seen. 
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Fig.  2  Plan- view  SEM  image  of  Au  NRs  of  sample  A  on  a  Si  substrate  after  liftoff. 
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Fig.  3  Normalized  extinction  spectra  (the  left  ordinate)  of  Au  NRs  measured  in  water  solution  of 
samples  A  and  B.  The  LSP  dipole  resonance  peaks  are  located  at  1240  and  1030  nm  for  samples  A  and 
B,  respectively.  The  light  source  spectrum  (the  right  ordinate)  of  the  used  OCT  system  is  also  shown. 
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Fig.  4  Normalized  extinction  spectra  of  Au  NRs  of  samples  A  and  B  obtained  from  simulations  when 
the  LSP  resonance  excitation  light  is  incident  along  the  z  axis  with  the  linear  polarization  in  the  y 
direction.  Two  sets  of  curves  are  shown  for  the  two  cases  of  1.33  (in  water)  and  1.4  (in  tissue)  in 
refractive  index  (numbers  shown  in  the  parentheses  of  the  labels).  The  light  source  spectrum  of  the 
used  OCT  system  is  also  shown.  The  simulated  NR  geometry  and  the  used  coordinate  system  are 

shown  in  the  insert. 


Fig.  5  Normalized  extinction  spectra  of  Au  NRs  of  samples  A  and  B  obtained  from  simulations  when 
the  LSP  resonance  excitation  light  is  incident  in  the  y  direction  with  the  linear  polarization  along  the  x 
(for  curves  A/ring  and  B/ring)  and  z  (for  curves  A/axis  and  B/axis)  axes. 
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Fig.  6  Simulation  results  of  spectral  dependencies  of  extinction  (Ext),  scattering  (Sea),  and  absorption 
(Abs)  cross  sections  of  a  single  Au  NR  for  samples  A  and  B. 
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Fig.  7  OCT  images  of  a  pig  adipose  sample  taken  before  (a),  0  min  (b),  30  min  (c),  and  60  min  (d) 
after  the  application  of  a  droplet  of  Au  NR  water  solution  of  sample  B.  The  sample  is  continuously 

scanned  with  the  OCT  system. 


after  the  application  of  a  droplet  of  Au  NR  water  solution  of  sample  A.  Part  (e)  shows  the  image 
similar  to  that  in  part  (d)  after  the  Au  NRs  accumulated  on  the  surface  is  washed  out.  The  sample  is 
continuously  scanned  with  the  OCT  system. 
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Fig.  9  OCT  images  of  (a)  a  pig  adipose  sample  with  delivered  resonant  Au  NRs  (a  duplicate  of  Fig. 
8(e))  and  (b)  a  fresh  pig  adipose  sample  (without  being  refrigerated  or  using  any  Au  NR). 


Fig.  10  OCT  scan  intensity  distributions  along  the  dashed  lines  in  Figs.  9(a)  (Adipose  with  NR)  and 

9(b)  (Fresh  Adipose). 
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Fig.  1 1  OCT  images  of  a  pig  adipose  sample  taken  before  (a),  0  min  (b),  30  min  ( 

c),  and  60  min  (d) 

after  the  application  of  a  droplet  of  Au  NR  water  solution  of  sample  A.  Part  (e)  shows  the  image 
similar  to  that  in  part  (d)  after  the  Au  NRs  accumulated  on  the  surface  is  washed  out.  The  sample  is 

scanned  once  every  30  min. 
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Depth  (jLim) 

Fig.  12  OCT  scan  intensity  distributions  along  the  vertical  lines  in  Figs.  11(a)  (Before)  and  11(e) 

(Wash  out). 
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Part  II:  Absorption  enhancement  of  an  amorphous  Si  solar  cell  through  surface 
plasmon-induced  scattering  with  metal  nanoparticles 

1.  Introduction 

In  an  amorphous  Si  (a- Si)  solar  cell,  because  of  the  low  carrier  mobility,  a  thinner  semiconductor  layer 
has  the  advantage  of  higher  carrier  harvest  in  generating  photocurrent.  Also,  the  degradation  issue  in 
an  a-Si  solar  cell,  known  as  the  Staebler-Wronski  effect  [1],  requires  the  minimization  of  the  active 
layer  thickness  to  avoid  carrier  loss.  However,  a  thin  semiconductor  absorbing  layer  may  result  in 
inefficient  absorption  and  low  device  efficiency.  Therefore,  it  is  important  to  use  certain  methods  for 
maintaining  efficient  absorption  under  the  condition  of  a  limited  semiconductor  thickness  such  that 
high  short-circuit  current  and  high  filling  factor  of  such  a  device  can  be  achieved.  Various  methods  for 
enhancing  the  absorption  of  a  thin  a-Si  solar  cell  have  been  proposed  [2-6].  Among  those  methods, 
sunlight  scattering  by  a  nanostructure,  which  is  deposited  on  the  top  surface  of  a  solar  cell,  has  shown 
to  be  an  effective  technique.  In  particular,  using  a  surface  metal  nanostructure  for  inducing  surface 
plasmon  (SP)  to  enhance  solar  cell  absorption  has  attracted  research  attention  in  various  photovoltaic 
materials,  including  organic  [7-15],  dye  [16,17],  GaAs  [18,19],  InGaN  [20],  and  Si  [21-30].  In  each  of 
those  implementations,  the  absorption  efficiency  in  a  certain  spectral  band  or  the  overall  solar  cell 
efficiency  was  more  or  less  enhanced.  The  possible  mechanisms  of  absorption  enhancement  with 
surface  metal  nanostructures  include:  1)  energy  transfer  from  excited  SPs  into  electrons  for 
contributing  to  photovoltaic  current  [9];  2)  increased  optical  path  of  incident  light  in  the  absorbing 
layer  through  SP  scattering  [19,30];  3)  increased  interaction  time  of  light  with  absorbing  material 
through  SP  (from  ~1  to  ~10  fs)  [21];  4)  enhanced  absorption  through  intense  SP  field  [20,21,30];  and 
5)  enhanced  light  incidence  into  the  solar  cell  material  through  the  forward  scattering  of  the  metal 
nanostructures  on  the  surface  [22,28].  Among  various  metal  nanostructures  for  the  aforementioned 
application,  the  use  of  metal  nanoparticles  (NPs)  for  inducing  localized  surface  plasmon  (LSP) 
resonance  represents  one  of  the  most  practical  approaches  because  of  fabrication  simplicity. 
Experimental  implementations  [19,23,25,26]  and  numerical  simulations  [20,23,26,29,30]  have 
demonstrated  the  effectiveness  and  practical  feasibility  of  this  approach.  However,  in  the  previous 
studies  of  the  solar  cells  with  Au  or  Ag  NPs  on  the  top  surfaces,  significant  reductions  of  sunlight 
absorption  were  observed  in  the  photon  energy  ranges  higher  than  the  corresponding  LSP  resonance 
energies,  leading  to  small  overall  increases  of  sunlight  absorption  or  solar  cell  efficiency  [18,23,29,30]. 
The  absorption  reduction  originates  from  the  sign  change  of  the  polarizability  of  the  used  metal  NP, 
which  may  cause  a  destructive  interference  between  the  incident  light  and  NP-scattered  light  in  the 
absorption  layer  of  the  solar  cell.  Although  the  use  of  dielectric  NP  on  the  device  top  surface  can 
avoid  the  destructive  interference  in  the  forward  direction,  the  overall  scattering  strength  is  weak  due 
to  the  small  polarizability  of  a  dielectric  NP  [28].  Therefore,  searching  an  appropriate  metal  type  for 
placing  on  the  solar  cell  top  surface  to  produce  constructive  interference  in  the  forward  direction  and 
hence  enhance  sunlight  absorption  is  an  important  issue  for  solar  cell  development.  In  choosing  metal 
type  for  such  application,  one  with  a  larger  LSP  resonance  energy  is  preferred  to  maximize  the 
forward  constructive  interference  effect  in  the  solar  spectrum  range.  Meanwhile,  since  the  geometry 
of  the  metal  NP  and  its  location  can  also  effectively  influence  the  scattering  behaviors,  these  factors 
also  deserve  detailed  investigation. 

In  this  part  of  report,  we  demonstrate  the  simulation  results  of  absorption  enhancement  in  an  a-Si 
solar  cell  by  depositing  metal  NPs  on  the  device  top  and  embedding  metal  NPs  in  a  layer  above  the  A1 
back-reflector  (back-contact).  From  the  evaluated  phase  distribution  results,  it  is  found  that  the 
absorption  increase  is  mainly  caused  by  the  enhanced  near-field  constructive  interference  effect  in  the 
forward  direction  through  the  mixture  of  LSP  resonance  and  Fabry-Perot  oscillation,  such  that  an 
increased  amount  of  sunlight  energy  is  distributed  in  the  a-Si  absorption  layer.  The  comparisons 
between  three  used  metals  of  Al,  Ag,  and  Au,  A1  NPs  show  the  most  efficient  absorption  enhancement, 
which  is  also  larger  than  that  of  using  Si02  NPs  of  the  same  geometry.  Also,  by  arranging  both  top  and 
bottom  metal  NPs  in  a  solar  cell,  further  absorption  enhancement  can  be  obtained.  With  a  vertically 
aligned  top  and  bottom  Al  nanocylinder  (NC)  structure,  the  solar  cell  absorption  can  be  increased  by 
52  %.  In  section  2  of  this  part  of  report,  the  problem  geometry  and  simulation  method  are  described. 
The  simulation  results  with  periodical  distributions  of  NPs  on  the  top  of  a  solar  cell  are  discussed  in 
section  3.  Then,  the  comparisons  between  the  cases  of  a  periodical  NP  distribution  and  a  single  NP  are 
presented  in  section  4.  Next,  the  effects  of  adding  a  bottom  metal  NP  distribution  are  reported  in 
section  5.  Finally,  the  conclusions  are  drawn  in  section  6. 
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2.  Problem  geometry  and  simulation  method 

Figures  1(a)  and  1(b)  show  the  solar  cell  structures  used  for  simulation  study.  Except  the  NPs,  the 
simulation  problem  geometries  represent  the  standard  structures  of  a-Si  solar  cells.  Either  device 
structure  consists  of  (from  the  bottom)  an  A1  layer  as  the  back  contact,  a  thin  n-type  a- Si  layer,  an 
intrinsic  a-Si  layer,  a  thin  p-type  a-Si  layer,  and  an  indium- tin- oxide  (ITO)  layer.  The  front  A1  contact 
is  neglected  in  simulation.  In  the  problem  geometry  of  Fig.  1(a),  which  will  be  referred  to  as  the  case 
of  “top  NP”,  a  metal  NP,  either  an  NC  or  a  nanosphere  (NS),  is  placed  at  the  center  of  the  top  surface 
in  the  shown  simulation  window.  The  simulation  window  stands  for  a  period  in  the  configuration  of  a 
periodical  metal  NP  distribution  in  both  the  x  and  y  directions.  To  demonstrate  the  LSP  effect,  the 
simulations  of  using  Si02  NC  and  NS  are  also  performed  for  comparison.  In  the  problem  geometry  of 
Fig.  1(b),  which  will  be  referred  to  as  the  case  of  “double  NPs”,  to  embed  the  bottom  metal  NP, 
another  ITO  layer  is  deposited  between  the  p-type  a-Si  layer  and  the  A1  back  contact.  The  NP  in  the 
bottom  ITO  layer  contacts  the  bottom  A1  layer.  To  show  the  location  of  the  bottom  metal  NP,  only 
one-half  the  simulation  window  is  shown  in  Fig.  1(b).  The  sunlight  is  assumed  to  be  normally  incident 
upon  the  device  from  the  top.  Although  the  sunlight  is  un-polarized,  without  loss  of  generality,  we  will 
consider  polarized  sunlight  with  the  polarization  oriented  in  the  x  direction.  When  the  NPs  are 
periodically  distributed  on  the  device  surface  and  such  a  periodical  arrangement  does  not  much  affect 
the  NP  scattering  behaviors,  or  when  the  NPs  are  randomly  distributed,  the  simulation  results  based  on 
polarized  sunlight  can  be  directly  applied  to  the  real  case  of  un-polarized  sunlight  due  to  the 
two-dimensionally  isotropic  nature  of  our  problem;  otherwise,  slight  modifications  are  needed.  In  Fig. 
1(a),  the  thicknesses  of  the  ITO,  p-type  a-Si,  intrinsic  a-Si,  and  n-type  a-Si  layers  are  assumed  to  be  20, 
10,  110,  and  10  nm,  respectively,  which  represent  a  widely  used  device  structure  of  this  kind  [29,30]. 
Here,  the  A1  back-contact  or  reflector  is  assumed  to  be  infinitely  thick.  In  Fig.  1(b),  the  thicknesses  of 
the  top  ITO,  p-type  a-Si,  intrinsic  a-Si,  n-type  a-Si,  and  bottom  ITO  layers  are  20,  10,  100,  10,  and  60 
nm,  respectively.  The  change  of  the  thickness  of  the  intrinsic  a-Si  layer  from  110  into  100  nm  means 
to  optimize  the  absorption  efficiency  of  the  solar  cell  when  NPs  are  not  applied  to  the  device.  The 
simulation  window  has  a  dimension  of  250  nm  in  either  x  or  y  direction.  In  the  simulations,  various 
metal  and  Si02  NCs  and  NSs  of  the  same  geometries  are  used.  The  NC  has  the  height  of  75  nm  and 
cross-section  diameter  of  100  nm.  The  NS  has  the  diameter  of  100  nm. 

The  simulation  computations  are  performed  with  the  finite-element  method  (FEM).  The 
perfectly  matched  layers  (PMLs)  are  used  in  the  +z  and  -z  directions  for  simulating  the  infinitely 
extended  air  and  A1  layers.  The  periodic  boundary  conditions  are  imposed  to  the  simulation  window  in 
the  +x,  -x,  +y,  and  -y  directions  to  model  the  problem  periodicity.  To  investigate  field  scattering  due  to 
the  NP,  we  use  the  well-known  scattered-field  formula  for  simulation,  in  which  the  total 
electromagnetic  field  is  divided  into  an  unperturbed  part  and  a  scattered  part  [31,32].  It  is  noted  that 
the  maximum  optical  intensity  in  the  absorbing  layer  of  the  solar  cell  can  be  obtained  when  the 
unperturbed  field  and  scattered  field  are  in  phase  in  the  forward  direction.  The  unperturbed  field  is  the 
electromagnetic  solution  with  the  designated  solar  cell  structure,  but  without  the  NP  under  the 
condition  of  an  incident  plane  wave.  It  can  be  obtained  from  an  explicit  analytic  expression.  The 
scattering  due  to  the  NP  is  then  described  by  the  scattered  field.  In  the  situation  with  NP,  the  total  field 
with  the  given  unperturbed  part  in  a  non-PML  region  is  governed  by  the  two  Maxwell’s  source- free 
curl  equations.  Because  a  PML  serves  to  absorb  the  scattered  field,  in  a  PML  region,  only  the  scattered 
part  needs  to  be  described.  Note  that  the  continuity  of  the  tangential  component  of  the  scattered  field 
is  required  at  the  interface  between  a  PML  region  and  a  non-PML  region.  Rigorously,  the  numerical 
treatment  in  the  FEM  is  implemented  by  using  the  corresponding  variational  description  [33].  In 
numerical  computations,  a  commercial  software  (COMSOL)  was  used.  It  has  provided  us  with  a 
reliable  framework  of  the  aforementioned  scattered-field  formulation.  The  accuracy  of  our  numerical 
results  was  verified  through  the  convergence  test.  The  minimum  division  mesh  size  was  2  nm  at  the 
metal/dielectric  interface.  The  maximum  mesh  size  of  an  individual  dielectric  region  is  set  to  be  one 
tenth  the  shortest  wavelength  in  material  of  our  concern  (300  nm  in  free  space).  The  dielectric 
constant  of  a-Si  is  obtained  from  the  measurement  data  in  [34]  except  that  the  absorption  coefficient 
near  the  band  edge  is  corrected  based  on  the  Forouhi-Bloomer  model  [35].  The  used  refractive  index 
and  absorption  coefficient  as  functions  of  wavelength  for  simulations  are  shown  in  Fig.  2.  The 
dielectric  constants  of  various  metals  are  obtained  from  [34].  Fixed  refractive  indices  of  1.46  and  1.7 
are  used  for  Si02  and  ITO,  respectively. 
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3.  Simulation  results  with  a  periodic  NP  distribution  on  the  top 

Figure  3  shows  the  photon  absorption  rates  as  functions  of  wavelength  with  Al,  Ag,  Au,  and  Si02  NCs 
on  the  top  of  solar  cells  (see  Fig.  1(a)).  For  comparison,  the  reference  case,  in  which  no  NP  is  used, 
and  the  photon  flux  of  AM  1.5G  are  also  demonstrated.  The  results  shown  in  Fig.  3  are  obtained  under 
the  condition  of  periodical  NP  distributions  with  the  period  of  250  nm  in  both  x  and  y  directions.  Here, 
one  can  see  the  oscillatory  behavior  of  the  reference  curve  due  to  the  Fabry-Perot  effect  in  the  vertical 
direction.  When  an  Au  or  Ag  NC  is  placed  on  the  top,  the  photon  absorption  in  the  UV-blue-green 
range  is  essentially  reduced.  However,  on  the  long-wavelength  side,  absorption  is  enhanced  except  in 
the  range  between  690  and  800  nm,  in  which  a  Fabry-Perot  peak  in  the  reference  case  becomes 
weakened  when  a  metal  NC  is  used.  As  a  result,  the  integrated  photon  absorption  rate  in  the  case  of 
Au  (Ag)  NC  is  reduced  (enhanced)  by  9  (10)  %,  when  compared  with  the  reference  case.  The 
integrated  photon  absorption  rates  under  various  conditions  are  listed  in  Table  I.  In  Fig.  3,  one  can  also 
see  that  when  an  Al  NC  is  placed  on  the  device  top,  the  photon  absorption  rate  is  always  enhanced  in 
the  whole  concerned  spectral  range  except  that  between  690  and  800  nm,  leading  to  the  increase  of 
integrated  photon  absorption  rate  by  39  %.  For  comparison,  a  Si02  NC  of  the  same  geometry  is  also 
used  for  evaluating  the  photon  absorption  rate.  In  this  situation,  the  photon  absorption  rate  is  also 
significantly  increased  (by  16  %),  as  compared  with  the  reference  case.  Because  the  absorption  layer 
of  the  solar  cell  is  below  the  NP  by  20-150  nm,  the  absorption  enhancement  is  supposed  to  be  due  to 
the  enhanced  near-field  constructive  interference  effect  of  sunlight  in  the  absorption  layer.  The 
interference  behavior  is  controlled  by  the  relative  phase  between  the  aforementioned  unperturbed  and 
scattered  fields.  The  phase  differences  between  the  unperturbed  and  scattered  fields  at  the  location 
right  below  the  NP  center  with  a  depth  of  65  nm  from  the  top  surface  of  the  p-type  a-Si  layer  as 
functions  of  wavelength  for  the  four  kinds  of  NP  are  shown  in  Fig.  4.  Here,  one  can  see  that  among 
the  metal  NPs,  the  phase  difference  of  the  Al  case  is  generally  smaller  than  those  of  the  Au  and  Ag 
cases.  However,  the  phase  difference  of  the  Si02  case  is  generally  the  smallest  among  the  four  cases. 
The  non-zero  phase  difference  in  the  Si02  case  is  mainly  due  to  the  layered  structure  of  the  solar  cell. 
The  deviations  of  the  phase  difference  curves  from  the  Si02  case  in  using  metal  NPs  are  strongly 
related  to  LSP  resonances.  According  to  the  Mie  theory  [36],  the  phase  distribution  of  the  perturbed 
field  behind  a  metal  NP  is  usually  significantly  altered.  However,  because  of  the  complex  solar  cell 
structure,  the  spectral  distributions  of  phase  difference  become  complicated  even  though  the 
extinction  spectral  peaks  of  the  Al,  Ag,  and  Au  NP  cases  can  be  identified  to  be  at  420,  570,  and  640 
nm,  respectively.  These  wavelengths  are  marked  by  the  vertical  arrows  in  Fig.  4.  It  is  noted  that  all  the 
four  curves  in  Fig.  4  have  peaks  around  660  nm.  Also,  the  three  metal  curves  have  a  common  dip 
around  720  nm,  which  is  shifted  to  -740  nm  in  the  Si02  curve.  The  peaks  and  dips  correspond  to 
Fabry-Perot  resonance  and  anti-resonance  spectral  positions,  respectively,  which  have  been  influenced 
by  the  NP.  Here,  the  common  Fabry-Perot  behaviors  lead  to  the  similar  spectral  features  beyond  630 
nm  in  wavelength.  It  is  noted  that  by  using  a  Si02  NP,  although  the  optical  intensity  in  the  forward 
direction  can  be  relatively  stronger  when  compared  with  those  in  other  directions,  its  overall  scattering 
is  not  as  strong  as  that  of  an  Al  NP,  which  has  the  LSP  resonance  effects,  such  that  its  absorption 
enhancement  is  smaller  than  that  of  the  Al  NP  case. 

The  results  in  Figs.  3  and  4  show  the  superiority  of  using  Al  NC  for  increasing  absorption 
through  the  enhancement  of  the  constructive  interference  effect  in  the  forward  direction.  To 
understand  the  effects  of  NP  geometry,  we  repeat  the  simulation  for  Fig.  3  by  replacing  NCs  by  NSs. 
The  results  similar  to  Fig.  3  are  shown  in  Fig.  5.  Here,  one  can  see  that  Si02  NS  also  significantly 
enhances  photon  absorption  rate  (also  by  16  %  when  compared  with  the  reference  case,  as  shown  in 
Table  I).  Al  NS  still  leads  to  the  largest  absorption  enhancement.  By  using  Al  (Ag)  NSs,  the  photon 
absorption  rate  can  be  increased  by  31  (22)  %,  when  compared  with  the  reference  case  (see  Table  I). 
Also,  Au  NS  results  in  slight  absorption  enhancement  (by  3  %).  In  the  Ag  and  Au  cases,  the  use  of  NS 
results  in  increased  enhancements;  however,  in  the  Al  case,  the  use  of  NS  leads  to  reduced 
enhancement,  when  compared  with  the  individual  cases  of  NC.  Figure  6  shows  the  phase  difference 
data  similar  to  those  in  Fig.  4  for  the  cases  of  NS.  Here,  one  can  see  the  similar  behaviors  to  those 
shown  in  Fig.  4.  The  Si02  case  has  generally  the  smallest  phase  difference  between  the  unperturbed 
and  scattered  light,  followed  by  the  Al  case.  The  three  wavelengths  for  the  maximum  extinctions  (LSP 
resonances)  at  310,  415,  and  550  nm  of  the  Al,  Ag,  and  Au  cases,  respectively,  are  also  marked  by 
vertical  arrows  in  Fig.  6.  The  significantly  shorter  wavelengths  for  LSP  resonances  in  the  cases  of  NS 
are  attributed  to  the  much  smaller  contact  area  of  an  NS  with  the  device  surface.  The  LSP  resonances 
in  the  NC  cases  are  strongly  influenced  by  the  contacting  ITO  layer.  It  is  noted  that  both  Si02  curves 
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in  Figs.  4  and  6  show  a  common  peak  near  480  nm,  around  which  the  curves  of  metal  NCs  and  NSs 
also  show  local  maxima.  This  feature  can  be  caused  by  the  interference  behavior  of  the  multi-layer 
solar  cell  structure  and  is  common  in  all  cases.  In  the  A1  case,  because  of  the  end-facet  contact  of  an 
NC  with  the  top  ITO  of  1.7  in  refractive  index,  the  resonance  wavelength  of  the  dominant  substrate 
LSP  mode  [32]  is  located  at  420  nm  (see  Fig.  4).  At  this  wavelength,  the  LSP-induced  absorption 
enhancement  is  more  prominent  due  to  stronger  sunlight  intensity,  when  compared  with  the  case  of  A1 
NS,  in  which  the  surrounding  air  results  in  the  LSP  resonance  wavelength  around  310  nm  (see  Fig.  6) 
of  weak  sunlight  intensity.  Therefore,  an  A1  NC  leads  to  a  higher  absorption  enhancement  than  an  A1 
NS.  On  the  other  hand,  in  the  cases  of  Ag  and  Au,  the  LSP  resonances  with  both  NC  and  NS  are 
located  in  the  spectral  range  of  high  sunlight  intensity.  However,  because  of  the  sharp  metal  edges, 
which  cause  higher  Ohmic  loss,  in  the  NC  geometry,  the  relatively  higher  dissipation  in  Ag  and  Au 
(compared  with  Al)  results  in  higher  loss  in  the  NC  geometry,  when  compared  with  the  NS  geometry. 
Hence,  the  absorption  enhancements  in  the  Ag  and  Au  NC  cases  become  lower  than  those  in  the 
corresponding  NS  cases. 

To  further  understand  the  enhanced  constructive  interference  effect  in  the  forward  direction,  the 
distributions  of  electrical  intensity  enhancement  ratios  with  respect  to  the  reference  case  for  the  total 
field  within  the  a-Si  regions  in  the  x-z  plane  are  shown  in  Fig.  7.  Here,  x  =  0  corresponds  to  the  center 
of  the  NP.  Also,  z  =  0  corresponds  to  the  bottom  of  the  n-type  a- Si  layer.  The  two  horizontal  white 
dashed  lines  represent  the  boundaries  between  the  p-type  and  intrinsic  a-Si  layers  and  between  the 
intrinsic  and  n-type  layers.  Figures  7(a)-7(h)  correspond  to  the  cases  of  Al  NC  at  525  nm  (enhanced), 
Ag  NC  at  600  nm  (enhanced),  Ag  NC  at  405  nm  (suppressed),  Si02  NC  at  525  nm  (enhanced),  Al  NS 
at  525  nm  (enhanced),  Ag  NS  at  525  nm  (enhanced),  Ag  NS  at  405  nm  (suppressed),  and  Si02  NS  at 
525  nm  (enhanced),  respectively.  Here,  one  can  see  from  Figs.  7(c)  and  7(g)  that  the  reduction  trend  of 
the  constructive  interference  effect  significantly  decreases  sunlight  absorption  of  the  whole  a-Si  layer. 
The  formation  of  the  two  bright  spots  in  the  situations  of  strong  constructive  interference  effect  is 
attributed  to  Fabry-Perot  oscillation.  It  is  noted  that  the  color  coding  scales  are  different  among 
different  cases  in  Fig.  7.  By  comparing  the  group  of  Figs.  7(a)  and  7(e)  with  the  group  of  Figs.  7(d) 
and  7(h),  one  can  see  that  the  scattering  strength  of  a  Si02  NP  is  indeed  lower  than  that  of  an  Al  NP, 
leading  to  the  relatively  smaller  absorption  enhancement  in  the  case  of  Si02  NP. 

4.  Simulation  results  with  a  single  NP  on  the  top 

The  results  given  above  correspond  to  the  condition  of  periodic  NP  distribution  on  the  device  top.  In 
this  section,  we  show  the  solar  cell  absorption  behaviors  when  a  single  metal  NP  is  used.  Figure  8 
shows  the  comparisons  of  photon  absorption  rate  as  a  function  of  wavelength  between  the  cases  of 
single  NC  and  periodic  NC  distributions  of  Al,  Ag,  and  Au.  Again,  the  curves  of  reference  and  AM 

I. 5G  are  also  depicted  for  comparison.  To  obtain  the  data  in  Fig.  8,  we  assume  that  the  effective 
absorption  region,  which  is  defined  for  evaluating  the  total  photon  absorption  rate,  is  a  circular  area  of 
282  nm  in  diameter.  The  assumption  of  this  diameter  value  means  to  make  the  effective  absorption 
area  the  same  as  that  (250  nm  x  250  nm)  in  the  cases  of  periodical  NP  distributions.  Here,  one  can  see 
that  the  spectral  oscillation  patterns  of  photon  absorption  rate  are  different  between  the  cases  of  single 
NC  and  periodic  NC  distribution.  The  periodic  distribution  of  NC  leads  to  certain  diffraction  effects 
and  changes  the  local  Fabry-Perot  behaviors.  Also,  the  LSP  interaction  between  neighboring  NPs  may 
also  cause  the  changes  of  spectral  patterns.  One  can  observe  alternative  variations  of  enhancement  and 
suppression  with  wavelength  in  comparing  the  cases  of  single  NC  and  periodic  NC  distribution.  The 
integrated  photon  absorption  rates  and  their  ratios  with  respect  to  the  reference  level  are  listed  in  Table 

II.  Here,  one  can  see  the  slight  increase  of  photon  absorption  rate  by  using  a  single  NP  of  any 
concerned  metal,  when  compared  with  the  condition  of  periodical  NP  distribution. 

Figures  9(a)  and  9(b)  show  the  distributions  of  electric  intensity  enhancement  ratios  over  that  of 
the  reference  case  at  500  nm  in  wavelength  in  the  x-y  plane  at  the  depth  of  65  nm  from  the  top  surface 
of  the  p-a-Si  layer  with  the  circles  centered  at  the  centers  an  Al  NC  (a)  and  an  Ag  NC  (b)  of  100  nm  in 
diameter  on  the  device  top.  The  values  of  D  in  nm  indicate  the  diameters  of  the  white  dashed  circular 
curves.  Here,  one  can  see  the  enhanced  and  reduced  intensity  distributions  right  below  the  Al  and  Ag 
NC,  respectively.  The  total  photon  absorption  rates  shown  in  Fig.  8  depend  on  the  definition  of  the 
effective  absorption  area.  Figure  10  shows  the  enhancement  ratios  of  photon  absorption  as  functions 
of  the  effective  absorption  diameter.  The  enhancement  ratio  is  defined  as  the  ratio  of  the  integrated 
photon  absorption  rate  in  the  case  with  NC  over  that  without  NC  (the  reference).  Therefore,  in  Fig.  10, 
when  the  effective  absorption  diameter  approaches  infinity,  the  enhancement  ratio  will  asymptotically 
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approach  one.  Here,  one  can  see  that  by  using  an  Au  NC,  when  the  effective  absorption  diameter  is 
between  200  and  500  nm,  the  enhancement  ratio  is  always  smaller  than  unity,  implying  that  in  this 
situation  the  incident  sunlight  is  either  absorbed  by  the  Au  NP,  backscattered  into  air,  or  scattered  in  a 
manner  that  the  photon  energy  is  pushed  away  from  the  designated  effective  absorption  region.  On  the 
other  hand,  by  using  an  A1  or  Ag  NC,  the  enhancement  ratio  decreases  monotonically  with  effective 
absorption  diameter.  In  particular,  the  use  of  an  Ag  NC  leads  to  a  faster  decay  into  unity,  implying  that 
the  effect  of  Ag  NC  scattering  is  limited  to  a  smaller  region  around  the  metal  NP  when  compared  with 
that  of  A1  NC  scattering.  By  using  an  A1  or  Ag  NC,  the  metal  NP  can  “attract”  sunlight  energy  into  the 
designated  effective  absorption  region  from  outside. 

5.  Comparison  between  the  cases  of  top,  bottom,  and  double  NPs 

The  results  above  show  the  enhanced  constructive  interference  effects  in  the  forward  direction  by 
placing  metal  and  dielectric  NPs  on  the  front  (top)  side  of  a  solar  cell.  In  this  section,  we  demonstrate 
the  effects  of  enhanced  backscattering  by  depositing  metal  NPs  on  the  backside  (bottom  side)  of  the 
solar  cell.  For  this  purpose,  we  consider  the  device  structure  shown  in  Fig.  1(b).  To  embed  a  metal  NP 
between  the  A1  back-contact  and  the  n-a-Si  layer,  an  ITO  layer  is  inserted  in  between.  Here,  we 
assume  that  an  A1  NC  is  embedded  in  this  ITO  layer.  For  our  simulations,  we  also  assume  that  the 
bottom  A1  NC  has  a  diameter  of  100  nm  and  a  height  of  60  nm  (smaller  than  that  of  the  top  metal  NC 
at  75  nm).  This  NC  height  is  the  same  as  the  thickness  of  the  bottom  ITO  layer  such  that  the  bottom 
metal  NC  contacts  with  both  the  n-a-Si  layer  above  and  the  A1  back-contact  below.  In  this  solar  cell, 
the  thickness  of  the  i-a-Si  layer  is  reduced  from  110  nm  to  100  nm  to  optimize  the  Fabry-Perot 
oscillation  for  maximizing  the  absorption  in  the  reference  case.  Figure  1 1  shows  the  photon  absorption 
rates  as  functions  of  wavelength  in  the  cases  of  periodical  top  A1  NC  distribution  (top),  periodical 
bottom  A1  NC  distribution  (bottom),  and  both  top  and  bottom  A1  NC  distributions  (double).  In  the  case 
of  double  NCs,  the  top  and  bottom  A1  NCs  are  vertically  aligned  (the  same  x-y  coordinates).  The 
curves  of  the  reference  case  (no  metal  NP)  and  AM  1.5G  are  also  shown  for  comparison.  Here,  one 
can  see  that  without  the  top  A1  NC,  the  use  of  the  bottom  A1  NC  suppresses  the  absorption  near  its 
peak  around  525  nm.  However,  two  absorption  peaks  are  formed  around  675  and  750  nm  by  using  the 
bottom  A1  NC  for  enhancing  the  integrated  absorption.  These  two  peaks  correspond  to  the  LSP 
resonances  of  the  bottom  A1  NC.  However,  the  formation  of  these  peaks  results  in  only  10  %  increase 
of  integrated  absorption,  which  is  significantly  lower  than  the  39  %  absorption  enhancement  in  the 
case  of  top  A1  NC  only,  as  shown  in  Table  III.  Nevertheless,  the  combination  of  the  vertically  aligned 
top  and  bottom  A1  NCs  can  lead  to  an  integrated  absorption  enhancement  of  52  %.  Compared  with  the 
case  of  top  A1  NC  only,  the  addition  of  the  bottom  A1  NC  results  in  a  significant  absorption  increase 
on  the  long- wavelength  side  even  though  that  around  the  peak  is  slightly  reduced.  In  this  situation,  the 
integrated  absorption  is  significantly  enhanced.  Figures  12(a)-12(i)  show  the  distributions  of  electrical 
intensity  enhancement  ratios  (over  that  of  the  reference  case)  in  the  x-z  plane  of  the  cases  in  Fig.  11 
with  (a)-(c)  for  the  cases  of  top  NC,  (d)-(f)  for  the  cases  of  bottom  NC,  and  (g)-(i)  for  the  cases  of 
double  NCs.  The  corresponding  wavelengths  are  525  nm  in  (a),  (d),  and  (g),  600  nm  in  (b),  (e),  and 
(h),  and  675  nm  in  (c),  (f),  and  (i).  The  two  horizontal  white  dashed  lines  represent  the  boundaries 
between  the  p-type  and  intrinsic  a-Si  layers  and  between  the  intrinsic  and  n-type  layers.  Here,  one  can 
see  that  at  each  wavelength,  the  intensity  enhancement  patterns  are  similar  among  the  cases  of  top  NC, 
bottom  NC,  and  double  NCs,  except  that  in  Fig.  12(c).  However,  the  levels  of  enhancement  are  quite 
different.  The  significantly  higher  enhancement  levels  in  the  case  of  double  NCs  can  be  clearly  seen. 

In  Figs.  11  and  12,  we  demonstrate  the  absorption  enhancement  when  the  top  and  bottom  A1  NCs 
are  vertically  aligned.  In  Figs.  13  and  14,  we  compare  the  results  of  different  misalignments  between 
the  top  and  bottom  A1  NCs.  Figure  13  shows  the  photon  absorption  rates  as  functions  of  wavelength  in 
various  cases  of  double  A1  NCs  with  the  vertical  alignment  shifted  horizontally  by  one-quarter  (A/4 
shift)  and  one-half  (A/2  shift)  the  period,  A,  in  the  x-  and  y  directions.  The  incident  sunlight  is 
assumed  to  be  x-polarized.  The  curves  of  no  shift,  reference,  and  AM  1.5G  are  also  shown  for 
comparison.  Here,  one  can  see  that  the  misalignments  cause  minor  modifications  of  the  absorption 
curves.  Clear  changes  are  observed  at  the  kinks  around  600  nm.  This  feature  can  be  due  to  the 
Fabry-Perot  effect  between  the  top  and  bottom  NCs  when  they  are  vertically  aligned.  As  the  two  NCs 
become  misaligned,  this  kink  becomes  less  prominent.  The  integrated  photon  absorption  rates  and 
their  ratios  with  respect  to  the  reference  level  of  various  cases  are  shown  in  Table  III.  Here,  one  can 
see  that  with  even  up  to  one-half  period  misalignments,  the  integrated  photon  absorption  rates  are  only 
slightly  reduced  (by  up  to  6  %).  Therefore,  the  vertical  alignment  of  the  top  and  bottom  NCs  is  not 
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critical  for  enhancing  absorption  under  the  assumed  condition  of  periodical  NC  distributions  at  the  top 
and  bottom.  Figures  14(a)- 14(1)  show  the  distributions  of  electrical  intensity  enhancement  ratios  (over 
that  of  the  reference  case)  in  the  x-z  plane  of  the  cases  in  Fig.  13  with  (a)-(c)  for  the  case  of  A/4 
shift-x,  (d)-(f)  for  the  case  of  A/2  shift-x,  (g)-(i)  for  the  case  of  A/4  shift-y,  and  (j)-(l)  for  the  case  of 
A/2  shift-y.  The  corresponding  wavelengths  are  525  nm  in  (a),  (d),  (g),  and  (j),  600  nm  in  (b),  (e),  (h), 
and  (k),  and  675  nm  in  (c),  (f),  (i),  and  (1).  The  two  horizontal  white  dashed  lines  represent  the 
boundaries  between  the  p-type  and  intrinsic  a-Si  layers  and  between  the  intrinsic  and  n-type  layers.  It 
is  noted  that  Figs.  14(g)- 14(1),  which  demonstrate  the  intensity  enhancement  ratio  distributions  in  the 
x-z  plane  by  shifting  the  alignments  in  the  y  direction,  look  similar  to  Figs.  12(g)-12(i).  Because  of  the 
symmetrical  NC  arrangements  between  the  x  and  y  axes,  the  misalignments  in  the  y  direction  does  not 
seem  to  affect  much  the  electrical  intensity  distributions.  However,  it  is  worth  mentioning  that  the 
slight  changes,  particularly  in  Fig.  14(j),  can  be  attributed  to  the  misalignment  in  the  other  direction 
and  the  oriented  incident  polarization  (in  the  x  direction).  Regarding  the  results  in  Figs.  14(a)- 14(f), 
we  can  first  notice  the  horizontal  shifts  of  the  bright  spots  in  Figs.  14(c)  and  14(f)  (at  675  nm  in 
wavelength)  from  the  center,  as  the  cases  shown  in  Figs  14(i)  and  14(1),  by  62.5  and  125  nm,  which 
correspond  to  A/4  and  A/2  shifts,  respectively.  As  mentioned  earlier,  the  feature  around  675  nm  in 
Figs.  11  and  13  is  due  to  the  LSP  resonance  of  the  bottom  A1  NC.  The  location  of  resonance  field 
distribution  shifts  with  the  NC  position.  The  bright  spots  in  Figs.  14(b)  and  14(e)  (at  600  nm  in 
wavelength)  also  show  horizontal  shifts  with  respect  to  those  in  Figs.  14(h)  and  14(k).  However,  the 
shift  ranges  are  significantly  smaller  than  the  corresponding  misalignment  ranges,  A/4  and  A/2, 
respectively.  As  pointed  out  earlier,  the  feature  around  600  nm  is  related  to  the  Fabry-Perot  effect 
between  the  top  and  bottom  A1  NCs.  Therefore,  although  they  are  separated  by  A/4  or  A/2,  the 
maximum  intensity  locations  (the  bright  spots)  in  between  do  not  shift  that  much  in  the  x  direction. 
Then,  at  525  nm  in  wavelength,  such  interactions  between  the  top  and  bottom  A1  NCs  are  stronger 
such  that  the  horizontal  shift  ranges  of  the  lower  bright  spots  become  even  smaller,  as  shown  in  Figs. 
14(a)  and  14(d),  when  compared  with  those  in  Figs.  14(g)  and  14(j).  The  detailed  interactions  between 
the  top  and  bottom  metal  NPs,  including  the  Fabry-Perot  effect  and  the  possible  LSP  coupling, 
deserve  further  investigation. 

6.  Conclusions 

In  summary,  we  have  demonstrated  the  simulation  results  of  absorption  enhancement  in  an  a-Si  solar 
cell  by  depositing  metal  NPs  on  the  device  top  and  embedding  metal  NPs  in  a  layer  above  the  A1 
back-reflector.  The  absorption  increase  resulted  from  the  near- field  constructive  interference  of  optical 
fields  in  the  forward  direction  such  that  an  increased  amount  of  sunlight  energy  was  distributed  in  the 
a-Si  absorption  layer.  Among  the  three  used  metals  of  Al,  Ag,  and  Au,  A1  NP  showed  the  most 
efficient  absorption  enhancement.  Between  the  two  used  NP  geometries,  Al  NCs  were  more  effective 
in  absorption  enhancement  than  Al  NSs.  Also,  a  random  distribution  of  isolated  NCs  could  lead  to 
higher  absorption  enhancement,  when  compared  with  the  cases  of  periodical  NC  distributions. 
Meanwhile,  the  fabrications  of  both  top  and  bottom  Al  NCs  in  a  solar  cell  resulted  in  further 
absorption  enhancement.  Misalignments  between  the  top  and  bottom  Al  NCs  did  not  significantly 
reduce  the  enhancement  percentage.  With  a  vertically  aligned  top  and  bottom  Al  NC  structure,  the 
solar  cell  absorption  could  be  increased  by  52  %. 
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Table  I  Integrated  photon  absorption  rates  and  their  ratios  with  respect  to  the  reference  level  of  various 
cases  in  Figs.  3  and  5.  The  photon  absorption  rate  of  AM  1 .5G  covers  the  whole  solar  spectral  range. 


AM 

1.5G 

Reference 

A1 

NC 

Ag 

NC 

Au 

NC 

Si02 

NC 

A1 

NS 

Ag 

NS 

Au 

NS 

Si02 

NS 

25.6* 

3.36  x  107 
s"1 

4.68* 

3.68* 

3.05* 

3.89* 

4.40* 

4.11* 

3.45* 

3.90* 

ratio 

1 

1.39 

1.10 

0.91 

1.16 

1.31 

1.22 

1.03 

1.16 

*  x  107  s'1 


Table  II  Integrated  photon  absorption  rates  and  their  ratios  with  respect  to  the  reference  level  of 
various  cases  in  Fig.  8.  The  photon  absorption  rate  of  AM  1.5G  covers  the  whole  solar  spectral  range. 


AM 

1.5G 

Reference 

Single 

A1 

Periodic 

A1 

Single 

Ag 

Periodic 

Ag 

Single 

Au 

Periodic 

Au 

25.6* 

3.36  x  107 
s'1 

4.76* 

4.68* 

3.84* 

3.68* 

3.10* 

3.05* 

ratio 

1 

1.42 

1.39 

1.14 

1.10 

0.92 

0.91 

*  x  107  s'1 


Table  III  Integrated  photon  absorption  rates  and  their  ratios  with  respect  to  the  reference  level  of 
various  cases  in  Figs.  11  and  12.  The  photon  absorption  rate  of  AM  1.5G  covers  the  whole  solar 
spectral  range. 


AM 

1.5G 

Reference 

Top 

NC 

Bottom 

NC 

Double  NCs 
or  no  shift 

A/4 

shift-x 

A/2 

shift-x 

A/4 

shift-y 

A/2 

shift-y 

25.6* 

3.38  x  107 
s'1 

4.68* 

3.71* 

5.13* 

5.04* 

4.97* 

5.04* 

4.93* 

ratio 

1 

1.39 

1.10 

1.52 

1.49 

1.47 

1.49 

1.46 

*  X  107  s'1 
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Fig.  1  (a)  Solar  cell  structure  in  a  simulation  window  with  an  NP  on  the  top;  (b)  Solar  cell  structure  in 
one-half  the  simulation  window  with  a  metal  NP  on  the  top  and  another  metal  NP  in  the  bottom  ITO 

layer  contacting  the  bottom  A1  layer. 
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Fig.  2  Used  refractive  index  and  absorption  coefficient  of  a-Si  for  simulations  obtained  from  Refs.  34 

and  35. 


Fig.  3  Photon  absorption  rates  as  functions  of  wavelength  with  Al,  Ag,  Au,  and  Si02  NCs  on  the  top  of 
a  solar  cell.  For  comparison,  the  reference  case,  in  which  no  NP  is  used,  and  the  photon  flux  of  AM 

1.5G  are  also  demonstrated. 
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Fig.  4  Phase  differences  as  functions  of  wavelength  between  the  unperturbed  and  scattered  fields 
when  Al,  Ag,  Au,  and  Si02  NCs  are  used  in  the  solar  cell  structure  shown  in  Fig.  1(a). 


Fig.  5  Photon  absorption  rates  as  functions  of  wavelength  with  Al,  Ag,  Au,  and  Si02  NSs  on  the  top  of 
solar  cells.  For  comparison,  the  reference  case,  in  which  no  NP  is  used,  and  the  photon  flux  of  AM 

1.5G  are  also  demonstrated. 


Fig.  6  Phase  differences  as  functions  of  wavelength  between  the  unperturbed  and  scattered  fields 
when  Al,  Ag,  Au,  and  Si02  NSs  are  used  in  the  solar  cell  structure  shown  in  Fig.  1(a). 
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Fig.  7  Distributions  of  electrical  intensity  enhancement  ratios  (over  that  of  the  reference  case)  within 
the  a-Si  regions  in  the  x-z  plane.  Here,  x  =  0  corresponds  to  the  center  of  the  NP.  The  two  horizontal 
white  dashed  lines  represent  the  boundaries  between  the  p-type  and  intrinsic  a-Si  layers  and  between 
the  intrinsic  and  n-type  layers.  Parts  (a)-(h)  correspond  to  the  cases  of  A1 NC  at  525  nm  (enhanced), 
Ag  NC  at  600  nm  (enhanced),  Ag  NC  at  405  nm  (suppressed),  Si02  NC  at  525  nm  (enhanced),  A1  NS 
at  525  nm  (enhanced),  Ag  NS  at  525  nm  (enhanced),  Ag  NS  at  405  nm  (suppressed),  and  Si02  NS  at 

525  nm  (enhanced),  respectively. 


Fig.  8  Comparisons  of  photon  absorption  rate  between  the  cases  of  periodical  top  metal  NC 
distributions  (curves  of  periodic  Al,  periodic  Ag,  and  periodic  Au)  and  single  top  metal  NP 
configurations  (curves  of  single  Al,  single  Ag,  and  single  Au).  The  curves  of  reference  (no  NP)  and 
photon  flux  of  AM  1 .5G  are  also  shown  for  comparison. 


Fig.  9  Distributions  of  electrical  intensity  enhancement  ratios  (over  that  of  the  reference  case)  at  500 
nm  in  wavelength  in  the  x-y  plane  at  the  depth  of  65  nm  from  the  top  surface  of  the  p-a-Si  layer  with 
the  circles  centered  at  the  centers  of  an  Al  NC  (a)  and  an  Ag  NC  (b)  of  100  nm  in  diameter  on  the 
device  top.  The  values  of  D  in  nm  indicate  the  diameters  of  the  white  dashed  circular  curves. 
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Effective  absorption  diameter  (nm) 


Fig.  10  Enhancement  ratios  as  functions  of  effective  absorption  diameter  with  respect  to  the  reference 
case  when  single  Al,  Ag,  and  Au  NCs  are  individually  placed  at  the  device  top,  corresponding  to  the 

data  in  Figs.  8  and  9. 


Fig.  11  Photon  absorption  rates  as  functions  of  wavelength  in  the  cases  of  periodical  top  Al  NC 
distribution  (top),  bottom  Al  NC  distribution  (bottom),  and  both  top  and  bottom  Al  NC  distributions 
(double).  In  the  case  of  double  NCs,  the  top  and  bottom  NCs  are  vertically  aligned.  The  curves  of  the 
reference  condition  (no  NP)  and  AM  1.5G  are  also  shown  for  comparison. 
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Fig.  12  Distributions  of  electrical  intensity  enhancement  ratios  (over  that  of  the  reference  case)  within 
the  a-Si  regions  in  the  x-z  plane  of  the  cases  in  Fig.  1 1  with  (a)-(c)  for  the  cases  of  top  NC,  (d)-(f)  for 
the  cases  of  bottom  NC,  and  (g)-(i)  for  the  cases  of  double  NCs.  The  corresponding  wavelengths  are 
525  nm  in  (a),  (d),  and  (g),  600  nm  in  (b),  (e),  and  (h),  and  675  nm  in  (c),  (f),  and  (i).  The  two 
horizontal  white  dashed  lines  represent  the  boundaries  between  the  p-type  and  intrinsic  a-Si  layers  and 

between  the  intrinsic  and  n-type  layers. 


Fig.  13  Photon  absorption  rates  as  functions  of  wavelength  in  various  cases  of  double  A1  NCs  with  the 
vertical  alignment  shifted  horizontally  by  one-quarter  (A/4  shift)  and  one-half  (A/2  shift)  the  period 
(A)  in  the  x-  and  y  directions.  The  incident  sunlight  is  assumed  to  be  x-polarized.  The  curves  of  no 
shift,  reference,  and  AM  1.5G  are  also  shown  for  comparison. 
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Fig.  14  Distributions  of  electrical  intensity  enhancement  ratios  (over  that  of  the  reference  case)  within 
the  a-Si  regions  in  the  x-z  plane  of  the  cases  in  Fig.  13  with  (a)-(c)  for  the  case  of  A/4  shift-x,  (d)-(f) 
for  the  case  of  A/2  shift-x,  (g)-(i)  for  the  case  of  A/4  shift-y,  and  (j)-(l)  for  the  case  of  A/2  shift-y.  The 
corresponding  wavelengths  are  525  nm  in  (a),  (d),  (g),  and  (j),  600  nm  in  (b),  (e),  (h),  and  (k),  and  675 
nm  in  (c),  (f),  (i),  and  (1).  The  two  horizontal  white  dashed  lines  represent  the  boundaries  between  the 
p-type  and  intrinsic  a-Si  layers  and  between  the  intrinsic  and  n-type  layers. 
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Part  III:  Surface  Plasmon  Coupling  with  a  Radiating  Dipole 


Surface  plasmon  (SP)  coupling  with  a  light  emitter  for  enhancing  light  emission  efficiency  is  an 
attractive  research  topic.  It  can  be  used  for  significantly  increasing  the  emission  strength  of  an 
orginally  weakly-emitted  molecule.  It  also  has  the  important  application  to  the  efficiency 
improvement  of  a  light-emitting  diode  (LED).  When  a  light  emitter  or  radiating  dipole  is  placed  near  a 
metallic  nanostructure,  the  induced  SP  resonance  can  couple  with  the  dipole  and  enhance  its  overall 
emission  efficiency.  It  is  supposed  that  such  a  coupling  process  can  increase  the  radiative  transition 
rate  of  the  radiating  dipole.  Based  on  the  SP  coupling  with  the  electron-hole  pairs  in  the  quantum 
wells  (QWs)  of  an  InGaN/GaN  QW  LED,  the  LED  emission  efficiency  improvement  has  been 
demonstrated  [1,2].  To  interpret  the  emission  efficiency  improvement,  numerical  studies  on  such  an 
SP-dipole  coupling  process  have  also  been  reported  [3,4].  In  such  studies,  various  numerical 
techniques  were  used  for  evaluating  the  enhanced  radiated  power  when  a  radiating  dipole  was  placed 
close  to  a  metal  surface  or  nanostructure.  However,  those  numerical  studies  were  based  on  the 
assumption  of  fixed  dipole  strength.  In  other  words,  the  efficiency  enhancement  is  simply  due  to  the 
additional  emission  of  the  equivalent  image  dipole,  represented  by  the  induced  SP  resonance.  Such  a 
coupling  picture  is  incomplete  because  the  induced  SP  resonance  can  influence  the  dipole  strength  or 
produce  a  feedback  effect  on  the  dipole  radiative  transition  rate.  A  more  complete  simulation  model  is 
needed  for  further  understanding  the  SP-dipole  coupling  scenarios.  In  this  letter,  we  numerically  study 
such  a  feedback  effect  by  assuming  a  two-level  system  for  the  dipole  oscillation  to  achieve  a  more 
complete  picture  of  SP-dipole  coupling.  It  is  found  that  the  dipole  emission  power  can  be  further 
enhanced  by  including  the  feedback  effect. 

In  this  study,  we  consider  a  system  consisting  of  a  metal  nanosphere  (NS)  and  a  nearby  radiating 
dipole,  as  shown  in  Fig.  1.  The  metal  NS  of  R  in  radius  and  <s^etai  in  dielectric  constant  is  centered  at 
the  coordinate  origin,  i.e.,  (0,0,0).  The  dielectric  constant  of  the  background  medium  is  £*out.  The 
radiating  dipole,  represented  by  p0(r,t)  =  p0(p0xx  +  p0yy  +  p0zz)e~i6,t  ^  is  located  at  (0,0, a).  The 

electric  field  of  the  radiating  dipole  induces  localized  surface  plasmon  resonance  (LSPR)  on  the 
metallic  NS,  resulting  in  a  strong  near  field  distribution  in  its  vicinity.  It  is  supposed  that  such  a  strong 
field  of  LSPR  can  affect  the  dipole  strength.  In  other  words,  the  LSPR  causes  a  feedback  effect. 

A  classical  dipole  corresponds  to  the  electron  transition  between  two  energy  levels.  Here,  we  use 
the  semi-classical  theory  for  solving  the  problem  of  electromagnetic  field  interaction  with  a  two-level 
system.  Since  the  size  of  the  NS  is  much  smaller  than  the  radiating  wavelength  and  the  interaction 
occurs  essentially  in  the  near  field  range,  the  quasistatic  approximation  is  used  and  the  induced 
electric  field  can  be  solved  analylically.  After  certain  derivations,  the  induced  electric  field  at  (0,0, a) 
can  be  written  as 

E (0,0, a)  =  pasE0fxx  +  p0yEJy  y  +  p0zEJz  z  (1) 

With 


En=- 


Po 


4£o*£0  Ka 


fX=fy=- 


(g  1  )f  Rs |  y  (»  + 1 )(»  +  2)  f  R  ^  , 


~  I 


s  +1  +  - 


1  u 


/.=(*“!)  T  I 


(n  +  2f 


£  +1  +  - 


n  + 1 
R 
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n  + 1 
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(2) 

(3) 

(4) 


g  =  gmetal(^)  (5) 

^out 

To  investigate  the  feedback  effect  of  the  LSPR  on  the  two-level  system,  we  solve  the  optical  Bloch 
equations  with  the  rotating  wave  approximation  and  the  dipole  approximation  [5].  For  mathematical 

simplicity,  we  define  _  =  P -  =PL>  -  =  P!_>  and  %  =  \pxf \fxf +\pv( \fv\2  +\pz\2 \fz\2 .  Then,  we 
Po  y  Po  1  Po 

can  derive  the  governing  equations  from  the  steady  state  solutions  of  the  optical  Bloch  equations  to 
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give 


\PoJxf  ,  \Poyfy\  ,  \PoJz\2 
\2  "v  Af  \  + r  as 


fl+  ] 

fl+  ^  ' 

1  ^  1 

l  1  +  B{) 

l  1  + 

1  ^  1 +  BS) 

and 


(6) 


(7) 


Pi  = 


Pm 


1  + 


4f, 


i+^(|a/J2+|a/v|2+Ia/z|2) 


(i=x,y,z) 


where  A  =  ^tL^X 
h  p0 


a  T  T 

,  B  =  -^i 


2^2  ?  Tj  is  the  spontaneous  decay  time,  T2  is  the  dephasing 


time,  and  p  is  the  dipole  matrix  element.  Once  Eq.  (6)  can  be  solved,  by  substituting  the  solution  into 
Eq.  (7),  one  can  obtain  the  resultant  dipole  components.  The  parameters  used  for  numerical 
calculations  include  77/=10‘9  s,  T2  =  2T7/,  p  =  5x1 0‘29  Cm,p0  =  2//,  and  eout=  6.25.  The  Drude  model  as 
co2 

^metai  (&>)  =  1 - - - is  used  for  the  calculation  of  frequency-dependent  dielectric  constant  of 

co(co  +  iy) 

the  metallic  NS.  In  the  following  calculations,  the  metal  is  assumed  to  be  silver  with  cop  =  1.19xl016 
rad/s  and  y  =  1.3 2x10 14  rad/s  and  the  radius  of  the  NS  is  set  at  10  nm. 

Two  orientations  of  the  dipole  are  considered,  including  one  along  the  z-axis  (the  z-dipole)  and 
the  other  along  the  x-axis  (the  x-dipole).  In  both  orientations,  four  locations  of  the  dipole,  i.e.,  a  =  45, 
50,  60,  and  70  nm,  are  considered  for  simulation.  The  dipole  strength  enhancement  ratios  as  functions 
of  wavelength  are  shown  in  Figs.  2  and  3  for  z-  and  x-dipole,  respectively.  Comparing  Figs.  2  and  3, 
one  can  see  that  the  enhancement  ratio,  enhanced  from  the  condition  of  no  metal  NS,  of  the  z-dipole  is 
larger  than  that  of  the  x-dipole.  The  enhancement  of  dipole  strength  is  about  35%  at  a  =  50  nm  for  the 
z-dipole  and  25%  at  a  =  45  nm  for  the  x-dipole.  Both  Figs.  2  and  3  clearly  show  two  peaks.  The  peak 
on  the  long  wavelength  side  is  caused  by  the  induced  dipole  resonance  on  the  metal  NS.  The  other 
peak  on  the  short  wavelength  side  corresponds  to  the  induced  quadrupole  resonance.  As  the  ratio  a/R 
becomes  larger,  the  shorter-wavelength  peak  diminishes  and  the  longer-wavelength  peak 
asymptotically  approaches  the  resonance  condition  of  a  metal  NS  excited  by  a  uniform  plane  wave 
[6]. 

The  total  radiated  power  of  the  dipole-NS  system  is  also  calculated  for  comparing  to  the  case 
without  the  feedback  effect.  Here,  the  total  radiated  power  is  normalized  by  the  corresponding  value 

2  f  3  4  \ 

of  an  isolated  dipole,  p  _  P o  .  n  ^lsp  ,  where  n  is  the  refractive  index  of  the  background 

l  3c3  ) 


medium,  which  is  set  at  2.5  by  assuming  that  the  dipole  and  NS  are  embedded  in  GaN,  and  colsp  = 
3.24xl015rad/s  is  the  FSPR  frequency  of  a  metal  NS  excited  by  a  uniform  plane  wave.  The  numerical 
results  are  shown  in  Figs.  4  and  5  for  the  z-  and  x-dipole,  respectively.  Here,  one  can  first  see  that  with 
the  feedback  effect,  the  total  radiated  power  at  either  dipole  or  quadrupole  resonance  is  significantly 
more  enhanced,  when  compared  with  the  case  of  no  feedback.  Also,  at  either  resonance,  the  total 
radiated  power  is  generally  more  enhanced  in  the  case  of  the  z-dipole,  when  compared  with  the  case 
of  the  x-dipole.  When  a  =  45  nm,  the  inclusion  of  the  feedback  effect  results  in  more  resonant 
enhancement  of  the  total  radiated  power  by  ~80  and  ~  72  %  for  the  z-  and  x-dipole,  respectively.  The 
total  radiated  power  enhancement  ratios  of  the  feedback  case  over  the  non-feedback  case  as  functions 
of  wavelength  with  the  z-  and  x-dipole  are  shown  in  Figs.  6  and  7,  respectively.  One  can  see  that  by 
including  the  feedback  effect,  the  enhancement  is  significantly  increased. 

In  summary,  we  have  demonstrated  the  simulation  results  of  the  enhanced  radiated  power  of  a 
radiating  dipole  when  it  was  coupled  with  the  FSPR  of  a  nearby  metal  NS  by  including  the  influence 
of  the  FSPR  on  the  strength  of  the  dipole,  which  was  represented  by  a  two-level  system  in  our 
semi-classical  model.  It  was  shown  that  with  this  feedback  influence,  the  total  radiated  power  could  be 
further  enhanced  by  more  than  70  %,  when  compared  with  the  condition  of  ignoring  the  feedback 
effect.  In  comparing  the  results  of  two  different  dipole  orientations,  it  was  discovered  that  at  a  fixed 
distance  from  the  metal  NS,  the  radial  dipole  orientation  led  to  more  significant  enhancement  of  total 
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radiated  power,  when  compared  to  the  orbital  dipole  orientation. 
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Fig.  1  configuration  of  the  simulation  system,  consisting  of  a  metal  NS  with  radius  R,  centered  at 
(0,0,0),  and  a  radiating  dipole  located  at  (0,0, a). 
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Fig.  2  Wavelength  dependencies  of  the  strength  of  the  z-dipole  at  various  dipole-NS  distances. 
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Fig.  3  Wavelength  dependencies  of  the  strength  of  the  x-dipole  at  various  dipole-NS  distances. 
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Fig.  4  Wavelength  dependencies  of  the  normalized  total  radiated  power  of  the  z-dipole  at  various 
dipole-NS  distances,  including  the  conditions  with  (F)  and  without  (NF)  the  feedback  effect. 
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Fig.  5  Wavelength  dependencies  of  the  normalized  total  radiated  power  of  the  x-dipole  at  various 
dipole-NS  distances,  including  the  conditions  with  (F)  and  without  (NF)  the  feedback  effect. 
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Fig.  6  Total  radiated  power  enhancement  ratios  of  the  feedback  case  over  the  non-feedback  case  as 
functions  of  wavelength  with  the  z-dipole. 
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Fig.  7  Total  radiated  power  enhancement  ratios  of  the  feedback  case  over  the  non-feedback  case  as 
functions  of  wavelength  with  the  x-dipole. 
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Future  Work: 


1.  Fabrication  of  metal  nanorings  of  various  geometries  and  understanding  their  localized 
surface  plasmon  behaviors 

We  will  continue  the  use  of  nanosphere  lithography  and  metal  secondary  sputtering  methods  for 
fabricating  Au  nanorings  of  various  ring  diameters  and  heights.  The  goal  is  to  control  the  localized 
surface  plasmon  resonance  wavelength  and  strength  for  covering  a  spectral  range  from  visible  through 
near-infrared.  Various  measurement  methods  will  be  used  for  understanding  their  resonance 
behaviors. 

2.  Comparison  of  localized  surface  plasmon  behaviors  between  Au  nanoring,  nanoshell,  and 
nanorod 

We  will  continue  the  cooperation  with  Professor  Randall  Lee  of  Houston  University  and 
Professor  C.  Y.  Mou  of  National  Taiwan  University  for  comparing  the  localized  surface  plasmon 
behaviors  of  our  Au  nanorings  with  their  nanoshells  and  nanorods.  The  effects  of  coated  porous  Si02 
layer  on  localized  surface  plasmon  resonance  behaviors  will  also  be  studied. 

3.  Experimental  study  of  surface  plasmon  coupling  with  a  dipole  through  cathodoluminescence 
measurement 

For  either  emission  or  absorption  (detection  or  imaging)  application,  the  fundamental  behavior  of 
surface  plasmon  coupling  with  an  emitting  or  absorber  needs  to  be  understood.  The  simplest  method 
for  this  study  is  photoluminescence.  However,  in  such  a  study,  the  excitation  strength  depends  on  the 
sample  structure  and  the  result  becomes  unreliable.  In  a  cathodoluminescence  measurement, 
luminescence  is  excited  by  incident  electron.  Such  electron  excitation  is  weakly  related  to  sample 
structure  for  more  accurately  understanding  the  coupling  behavior. 

4.  Further  numerical  study  of  surface  plasmon  coupling  with  a  dipole 

We  will  study  the  coupling  behaviors  of  one  surface  plasmon  mode  with  more  than  one 
dipole  to  see  the  effects  of  coherent  coupling  and  emission.  If  such  a  multiple-dipole  cohere 
nt  coupling  can  be  implemented,  even  stronger  enhancements  of  emission  and  absorption  throu 
gh  surface  plasmon  coupling  can  be  implemented. 

5.  Further  numerical  study  of  surface  plasmon  enhanced  solar  cell 

We  will  further  numerically  study  the  absorption  enhancement  effect  of  a  silicon  solar  cell  with 
metal  nanostructures  to  optimize  the  design  of  such  a  solar  cell.  Hopefully,  we  can  receive  equipment 
for  implementing  such  a  silicon  solar  cell  in  the  near  future. 
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